Superconducting pairing of spin polarons in the t-J model 
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The pairing of quasiparticles in a CuC>2 plane is studied within a spin polaron formulation of 
the t-t -J model. Our numerical solution of the Eliashberg equations unambiguously shows d-wave 
pairing between spin polarons on different sublattices mediated by the exchange of spin-fluctuations, 
and a strong doping dependence of the quasiparticle bandwidth. The transition temperature T c is 
an increasing function of J/t and crosses a maximum at an optimal doping concentration 8 op t. For 
the t-J model with J/t = 0.4 we obtain T c ~ 0.013* at S opt ~ 0.2. 
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Recent experimental evidence in favor of a ci-wave su- 
perconducting pairing in high-T c cuprates [|J supports 
theoretical studies of models with strong electron corre- 
lations @ . The minimal model describing hole motion in 
CuC>2 plane is the t-J || or t-t' -J H model. Numerical 
studies jHH^| of small t-J clusters suggest a d-wave su- 
perconducting instability. Yet to elucidate the nature of 
this pairing, an analytical treatment of the t — J model is 
needed. For this purpose we use a spin polaron formula- 
tion for the t—J model deduced in the region of small hole 
concentrations. A number of studies of this model [QJ|] 
predict that doped holes dressed by strong antiferromag- 
netic spin fluctuations propagate coherently as quasipar- 
ticles (spin polarons) with weight Z k — J/t. ft is quite 
natural to expect that the same spin fluctuations induce 
superconducting pairing of the spin polarons. Recently 
this problem has been treated in the framework of the 
standard BCS formalism assuming a rigid band model 
for the quasiparticles However, since the spin- 

fluctuation energy is of the same order as a quasiparticle 
bandwidth ~ J a strong coupling approach is necessary. 

In this paper we present the first consistent solution of 
the strong coupling spin polaron model at finite temper- 
atures and hole concentrations for normal and supercon- 
ducting states. A numerical solution of a self-consistent 
system of equations for hole and magnon Green functions 
proves singlet d-wave superconducting pairing. The gap 
function shows interesting additional structure on top of 
the simple A k = Ac^cos/cz — cosk y ) which reflects the 
Fermi surface geometry. The doping dependence of T c 
around S op t has the form of an inverted parabola, similar 
to experiment, and a T™ ax ~ Combining these re- 

sults with already existing weak coupling studies for the 
Hubbard model we argue that the spin-exchange 

pairing is the true mechanism for high-temperature su- 
perconductivity as proposed earlier by several groups 
based on more phenomenological approaches 
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antiferromagnetic (AF) background which has been suc- 
cessfully tested in the single hole case (tJ, ||. Spinless 
fermion operators hf and ft are introduced for holes 
on different sublattices, i.e. on the f (j)-sublattice the 
constrained electron operators Ci a = c,* CT (l — n%-tr) of 
the t-J model arc replaced by cijj = hf, = hf Sv" 
( s a = ft , c iT = ftSr), where Sf = Sf ± Sf are spin 
operators. This representation excludes doubly occupied 
states and takes into account strong AF spin correlations 
in the electron hopping. 

By employing the linear spin-wave approximation 
in terms of the Holstein-Primakoff operators: — 
di, (i €t), St — bf , (i El) and performing the 
Bogoliubov canonical transformation: ak = Vka/c + 



model: 



u k a_ 



ki 



we obtain the spin polaron 



Ht-.J = ^204 ~fk-q[g(k, q)a q + g(q - k, q)f3± q ] + h.c. 



kq 



zZ^(hth k + f+f k ) + Y J "M a « + ftp*)- (!) 



We will study a spin polaron model on a two sublattice 



Here g(k, q) = [ztj y/N/2)(u q j k - q + v q j k ) is the holc- 
magnon interaction, z — 4 is the number of the near- 
est neighbors on a square lattice with N sites, u k — 
{{l + Vk )/2 Vk ) 1/2 , v k = -sign(7 fc )((l-^)/2^) 1/2 , 
v k = yj 1 — 7j? , j k — ^(cosfej. + cosky). The next near- 
est neighbor hopping energy is e k = {At' cos k x cos k y — 
fi). The chemical potential /j, should be calculated self- 
consistently as a function of a hole concentration d and 
temperature T from the equation: S = (hfhi) + (ftfi)- 
The spin-wave energy is u> q = SzJ(l—5) 2 v q where (1— S) 2 
is the mean field renormalization factor. We neglect here 
the contact hole-hole interaction which is unimportant 
in the polaron pairing fic|| . The summation over wave- 
vectors in (Q) and below is restricted to N/2 points in 
the AF Brillouin zone. 

To discuss singlet superconducting pairing within the 
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spin polaron model (|lj), we consider the matrix Green 
function (GF) for holes on two sub-lattices G n h(k,z) = 
{(K I h k )) z = ((f£ | f k )) z and the anomalous GF 
G hf (k,z) = ((h+ | f+ k )) z = -((f+ k | h+)) z , where 
Zubarev's notation fl7j for the anticommutator GF was 
used with z = lu + ie. To obtain self-consistent equa- 
tions for these GF's we employ the self-consistent Born 
approximation (SCBA) which provided good results for 
the one-hole spectrum in the normal state In 
SCBA we get for the self-energies 



T, hh (k,ico n ) = -T^G hh (q,iuj m )Xl 1 j,_ q (uj n -uj m ), (2) 

q , rn 

T, h f(k,ioj„) = -T^G hf {q,iuj m )X 1 k 2 k ^ q (uj n -u m ). (3) 



where the Matsubara frequencies uj n — nT(2n + 1). The 
interaction functions are 
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FIG. 1. The density of states (DOS) for the hole concen- 
tration 5 = 0.06 at T = 0.012*. The solid (dashed) line cor- 
responds to calculations with the full (zero order) magnon 
Green function. In the inset DOS is given for S — 0.1, 0.25, 
0.35 (from right to left) for the zero order magnon spectra. 



*%q(vv) = 9 2 (k, q)D(q, -iu> v ) + g 2 (q - k, q)D(-q, iu) v ), 
>?k,q(uv) = g(k, q)g{q - k, q){D(q, -iiO v ) + D(-q, iu;^)}. 

The diagonal magnon GF D(q,u>) — ((a q \ af)) w in 
the zero order approximation is given by D°(q,uj) — 
(uj — Wq) -1 with the doping dependent magnon energy 
Lo q . The full magnon GF is determined by the matrix 
equation l) _1 ((7,aj) = (D°)~ 1 (q, lu) — tl(q,tu) where the 
renormalization of the magnon energy due to particle- 
hole excitations is described by the polarization operator 
n(<Z, uj). This is calculated in one- loop approximation 
using the fully renormalized hole-GF. 

The superconducting temperature T c is calculated 
from the linearized form of the Eliashberg equation for 
the gap-function 



4>(k,iu n ) = J^J^ ^l 2 k- P ( iuj n - iuJ m )G hh (p,iLj m ) 

p rn 

xGhh(-p,-ium)(t>(p,iv m ). (4) 



The first step is a self-consistent calculation of the nor- 
mal GF G h h(k,iu! n ) = (iu n + e k - S /lh (fc, icj„)) _1 with 
the self-energy operator (2) for a given concentration of 
1 2T 

holes 5= - +-^EfcEn G WM^n)- 



The numerical calculations were performed using fast 
Fourier transformation (FFT) Jl8| for a mesh of 64 x 64 k- 
points in the full Brillouin zone (0 < k Xl k y < 1), in units 
of 2-7T. In the summation over the Matsubara frequen- 
cies we used up to 200-700 points with constant cut-off 
^max = 10t. The FFT for the momentum integration is 
possible due to the particular momentum dependence of 
g(k, q). Usually 10 - 30 iterations were needed to obtain 
a solution for the self energy with an accuracy of order 
0.001. Pade approximation was used to calculate the hole 

spectral function A(k,ui) = Im ((h k | ht))u)+ie and 

the density of states (DOS) A(ui) on the real frequency 
axis. In Fig. 1 results for A{lo) of the t-t'-J model are 
shown for various doping concentrations. The peak in 
the DOS of width AW < J near the chemical poten- 
tial /x = results from the shallow quasiparticle disper- 
sion E(k) along the AF-zone boundary (Fig. 2(a)). We 
find that the shape of the quasiparticle dispersion even 
at S ~ 0.25 is still similar to the shape of the disper- 
sion in the single hole case. Yet a rigid band description 
fails since the scale AW and the total quasiparticle band- 
width W grow significantly with S . The peak of the DOS 
coincides with \x at the crossover from hole to electron 
like Fermi surfaces (FS). This occurs at a characteris- 
tic concentration which depends on t' . In Fig. 2(b) the 
FS at S = 0.25 is shown for the two models studied in 
this paper, the t-J and the t-t'-J model with t 1 = — O.lt. 
Our unit of energy is t = 1 (in reality t ~ OAeV for 
CuO-2 planes) and J/t = 0.4. The crossover from hole- 
to electron-like FS is consistent with the variation of the 
Hall constant in La x Sr2~ x CuO^. 
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FIG. 2. (a) The quasiparticle spectrum E(k) and (b) the 
Fermi surface (FS) E(k F ) = of the t - t - J{t - J) model 
for 8 — 0.25 is given by the solid (dashed) line. 

The particle-hole renormalization of the magnon prop- 
agator D leads to an instability at small q indicating the 
disappearance of AF long range order. In Fig.l we com- 
pare the DOS A(ui) at 5 = 0.06 calculated with D and 
Dq, which is not ill behaved. The small-g instability has 
only small effects on A(k, uj) and A(oj) since in the small 
q-regime the spin-charge coupling is small. Therefore we 
performed our calculations at higher 5 with Dq. Our 
main assumption here is, that the spin polaron approach 
gives a reliable description also in the spin liquid regime 
provided the AF correlation length is sufficiently large 
compared to the Cooper pair and polaron radius. The 
latter quantity is 2 lattice constants for J/t = 0.4 jl9| . 

The momentum dependence of the gap function 
A(k,uj = 0), A(k,u) = <fi(k, uj)/Z(k, uj), is shown in 
Fig. 3(a) for S = 0.25 and T/T c « 0.8. Here Z(k,uj) 
is an analytical continuation of the Eliashberg function 
Z(k,iuj n ) = (1 — IrriT,(k,iuj n ) / uJn)^ 1 . The gap function 
has the typical cZ-wave symmetry with two ridges result- 
ing from sharp changes of the interaction function at the 
FS. In Fig. 3(b) the frequency dependence of i?eA(k, uj) 
is shown for a set of (k x ,k y ) points marked in Fig. 2b: 
(1) inside the FS, (0, 0.19), (2) at the AF-zone bound- 
ary, (0.31, 0.19), (3) near the FS, (0.38, 0.19). The gap 



function changes sign after crossing the k x — k y = 0.19 
point where it is equal to zero. It is interesting that the 
characteristic energy cutoff for the pairing theory, which 
is of order J ~ 0.4 away from the FS (curve 1 ), be- 
comes much smaller near the FS (curves 2 and 3). The 
sharp change of the real part and the quite large values 
of imA(fc, uj) near the FS differ from the results for con- 
ventional superconductors. Since the Fermi energy Ep 
is of the order of the exchange energy J all quasiparti- 
cles contribute to the pairing state contrary to the weak 
coupling case in conventional superconductors. 




FIG. 3. (a) The gap function A(k, ui — 0) versus k (in units 
of 2n/a) and (b) Re A(k, ui) (Im A(k,tj) in the inset) versus 
ui for a set of (k x , k y ) points shown in Fig. 2b for t-t -J model 
( S = 0.25 and T/T c « 0.8.) 

The transition temperature T c is determined as the 
temperature where the highest eigenvalue of the lin- 
earized Eliashberg equation becomes unity. In all cases 
the symmetry of the corresponding eigenfunction (f>(k,uj) 
is d x 2_ y 2. In Fig. 4 the dependence of T c on hole con- 
centration is shown for t 1 = — O.li and t' — 0. These 
results are quite different from the monotonic increas- 
ing T c obtained within the weak coupling limit of the 
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BCS equation in |10| and the maximum of T c found in 
near half filling. In our case the maximum of T c at 
S ~ 0.25 (or at S ~ 0.20 for t' — 0) results from the Fermi 
level crossing of the peak in the density of states which 
coincides with the change of the FS topology. 
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FIG. 4. The superconducting temperature T c versus hole 
concentration 8 for J = 0.4, t' = — O.lt (solid line) and t' = 
(dashed line). 



culations are based on a two sublattice representation, 
which is suggested to provide a reasonable description of 
spin polaron quasiparticles and their pairing even in the 
spin liquid regime, i.e. for hole concentrations where the 
AF correlation length is sufficiently large. 

An important difference between the phcnomcnologi- 
cal spin-fluctuation theory and our approach is that pair- 
ing is dominated in the former by q ~ (jr, 7r) scattering 
and energy transfers AE < 50mev, whereas in our cal- 
culations high-energy spin-fluctuations with q near the 
AF-zone boundary are most important. Higher energy 
neutron scattering data in this momentum and energy 
range is however not yet available. 
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We have also studied the dependence of T c on the ex- 
change energy for J < 4. T c increases with J and sat- 
urates at T c ~ 0.025t for J ~ 3. However, we have not 
obtained a large drop of T c for J > 3 observed in small 
cluster calculations near phase separation [0, which is 
beyond the scope of our study. 

In summary, we have solved numerically Eliashberg 
equations for the strong coupling spin-polaron model. 
We have calculated the quasiparticle spectrum of spin 
polarons in the normal state and shown that they un- 
dergo superconducting d-wave pairing mediated by spin 
fluctuations. The high values of superconducting temper- 
ature and its doping dependence T c {5) is explained by a 
large peak in the density of states of the spin polaron 
quasiparticles in the vicinity of the chemical potential. 
A key difference from the van Hove scenario, however, is 
that the quasiparticle density of states has a width given 
by the interaction energy J which is similar to the pair- 
ing energy. We have found unconventional behavior for 
the d-wave gap function (a sharp change with energy and 
large damping near the FS) which suggests an explana- 
tion for some of anomalous properties of cuprate super- 
conductors observed in tunneling experiments (v-shape 
gap and large imaginary part), infrared absorption (no 
visible gap or gapless superconductivity) , ARPES (a line 
of gap nodes along (7r,7r) direction etc. Our cal- 
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